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ABSTRACT 



Recent models of spectral formation in magnetars called renewed attention on 
electron-photon scattering in the presence of ultra-strong magnetic fields. Investiga- 
O ' tions presented so far mainly focussed on mildly relativistic particles and magnetic 

> . scattering was treated in the non-relativistic (Thomson) limit. This allows for consis- 

' tent spectral calculations up to a few tens of keVs, but becomes inadequate in mod- 

elling the hard tails (< 200 keV) detected by INTEGRAL from magnetar sources. In 
this paper, the second in a series devoted to model the X-/soft 7-ray persistent spec- 
^ . trum of magnetar candidates, we present explicit, relatively simple expressions for the 

' magnetic Compton cross-section at resonance which account for Landau-Raman scat- 

tering up to the second Landau level. No assumption is made on the magnetic field 
strength. We find that sensible departures from the Thomson regime can be already 
. present at B 5 x 10^^ G. The form of the magnetic cross section we derived can be 

' easily implemented in Monte Carlo transfer codes and a direct application to magnetar 

. spectral calculations will be presented in a forthcoming study. 

' Key words: Radiation mechanisms: non-thermal - stars: neutron - X-rays: stars. 
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1 INTRODUCTION 



The recent discovery with the INTEGRAL satellite of hard X-ray tails in the (persistent) spectra of the magnetar candidates 
(the anomalous X-ray pulsars, AXPs, and the soft 7-repeaters, SGRs; e.g. Moroghctti 2002,) provides evidence that a sizeable 
fraction (up to ^ 50%) of the energy output of thes e sources is emitted above ~ 20 keV. Up to now, high-energy emission has 
been detected in two SGRs, 1806-20 a nd 1900-1-14 jMereghet ti ct al."2005': 'Cotz et al."2006''l. and three AXPs, 4U 0142-f614, 
IRXS J170849-4009 and IE 1841-045 jKu i per. Hermsen fc M cndcz 2004; Kuipcr ct al. 2006, see also for an updated summary 
of INTEGRAL observations ICotzlboOsi ) 14 This result come somehow unexpected, since the spectra of SGRs/ AXPs in the 



soft X-ray range (~ 0.1-10 keV) are well described by a two component model, a blackbody at kT ~ 0.2-0.6 keV, and a 
rather steep power-law with photon index Tgoft ~ 1.5-4. INTEGRAL observations have shown that in SGRs the power-law 
tail extends unbroken (or steepens) in the ~ 20-200 keV range, Thard ^ ^aoft ~ 1-5. In AXPs, which have steeper spectra in 
the soft X-ray band, a spectral upturn appears, i.e. the high-energy power-law is harder than the soft one, Vhard ~ 1 while 



Pso/t ~ 3-4 (see e.g. iRea et al.l 2008, for a joint spectral analysis of XMM-Newton and INTEGRAL data) 



Within the magnetar scenario, the persist ent 0.1-10 keV emission of SGRs and A XPs has been successfully interpreted 



in terms of the twisted magnetosphere model IjThompson. Lvutikov fc Kulkarnill2002l ). In an ultra-magnetized neutron star 



the huge toroidal field stored in the interior produces a deformation of the crust. The displacement of the footprints of the 
external (initially dipolar) field drives the growth of a toroidal component which, in turn, requires supporting currents. Charges 



* E-mail: nobili @pd.infn.it (LN); turolla . @pd.in fn.it (RT); sz@mssl.ucl.ac.uk (SZ) 

1 Very recently [Eevder. Walter fc Rauwl l|2008l ) reported the INTEGRAL detection of the AXP IE 1048.1-5937, but no spectral infor- 
mation are available yet. 
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flowing in the magnetosphere provide a large optical depth to resonant cyclotron scattering (RCS) and repeated scatterings 
of thermal photons emitted by the star s urface may then lead to the formation of a pow e r-law t ail. The original picture 



by Thompson. Lvutikov fc Kulkarni 
( 2007f l and iNobili. TuroUa fc Zan3 ' 



(I2OO2I) has been further explored b vlLvutikoy fc Gavriill (|2006l ). Fernandez fc ThompsonI 



2008. hereafter paper I). Recently Rea et al.l ( 2008t ) presented a systematic application 



Lvutikov fc Gavriill (|2006l l to a large sample of magnetars X-ray spectra finding a good 



of the ID, analytical RCS model of 
agreement with data in the 0.1-10 keV range. As shown in paper I, more sophisticated 3D Monte Carlo calculations of RCS 
spectra successfully reproduce soft X-ray observations. 

Despite the twisted magnetosphere scenario appears quite promising in explaining the magnetars soft X-ray emission, 
if and how it can account also for the hard tails detected with INTEGRAL has not been unambiguously shown as yet. 
Thompson fc Belobodor"wl l|2005l ) suggested that the hard X-rays may be produced either by thermal bremsstrahlung in 
the surface layers heated by returning currents, o r by synchrotron emission from pairs created higher up (~ 100 km) in 
the magnetosphere. Baring fc Hardine 1 2007 , 20081 ) have recently proposed a further possibility, according to which the soft 
7-rays may or iginate from resonant up - scattering of seed photons on a population of highly relativistic electrons. Previous 
investigations ( Lvutikov fc GavriiJbood : Fernandez fc ThompsonlboOTi : paper I) mainly focussed on mildly relativistic particles 
and magnetic scattering was treated in the non-relativistic (Thomson) limit. This is perfectly adequate in assessing the spectral 
properties up to energies <^ nieC'^ /j (here 7 is the typical electron Lorentz factor) since i) the energy of primary photons is 
low enough (~ 1 keV) to make resonant scattering onto electrons possible only where the magnetic field has dropped well 
below the QED critical field, ii) up-scattering onto electrons with 7 ~ 1 hardly boosts the photon energy in the hundred 
keV range, so electron recoil is not important. However, some photons do actually gain quite a large energy (because they 
scatter many times on the most energetic electrons) and fill a tail at energies > 50 keV. We caveat that, despite in previous 
works spectra have been computed up to the MeV range, they become untrustworthy above some tens of keV and can not 
be used to assess the spectral features that can arise due to electron recoil effects (i.e. a high-energy spectral break). Proper 
investigation of the latter demands a complete QED treatment of magnetic Compton scattering. This is mandatory if highly 
relativistic particles are considered because a photon can be boosted to quite large energies in a single scattering and, if it 
propagates towards the star, it may scatter again where the field is above the QED limit. 

Monte Carlo numerical codes for radiation transport in a magnetized, scattering medium, as the one we presented in 
paper I, make an excellent tool to investigate in detail the properties of RCS in the case in which energetic electrons are 
present in addition to the mildly relativistic particles which are responsible for the formation of the soft X-ray spectrum. 
The Co mpton cross-section for electron sc a tterin g in the presence of a magnetic field was first studied in the non rel ativistic 
limit by Canuto. Lodengu ai fc Rudcrman (Il971), and the QED expressi on was derived long ago b y many authors I Herold 
I979I : Daughertv fc Harding 1986.. : Bussard. Alexander fc Meszaroa 1986, : .Harding fc Daughertvlll99ll ). However, its form is so 



complicated to be often of little practical use in numerical calculations. Moreover, because of their inherent complexity, many 
of the published expressions are affected by misprints and the comparison between different derivations is often problematic. 
On the other hand, the use of the full expression of the cross section is especially needed when a detailed model of cyclotron line 
form ation has to be evaluated, including expectations for the line profile (see e.g. Daughertv fc Venturall9'73 : Arava fc Hardin3 
I999I ). In the situation we are considering, it is reasonable to assume that nearly all photons scatter at resonance. Non-resonant 



scattering contributions have negligible effects on shaping the overall spectrum, except possibly in the very neighborhood of 
a cyclotron line peak. 

Motivated by this, we present here explicit, relatively simple expressions for the magnetic Compton cross-section at 
resonance that can be then included in Monte Carlo calculations such as that described in Paper I. In doing so, we investigate 
the behaviour of the different terms and assess their relative importance. The final expressions have been cross-checked by 
comparing different published formulations, when specialized at resonance. A direct application of the results discussed here to 
spectral calculations will be presented in a forthcoming study (Nobili, Turolla fc Zane, in preparation). The paper is organized 
as follows. In ^we formulate the problem and summarize the main ingredients needed for a Monte Carlo simulation. In ^ 
we compute the relevant cross sections, specified at the resonance, while the transition rates (which are related to the natural 
width of the excited resonant levels) are given in |4l The creation of photon via spawning effects is discussed in ^ while 
36] contains a brief comparison between scattering and absorption. In ^ we compute the optical depth, which yields the 
probability of scattering. Conclusions follow in |8l 



2 FORMULATION OF THE PROBLEM 

In this section we briefiy outline our approach, discuss our working assumptions and introduce the mathematical notation. In 
dealing with relativistic scattering it is convenient to express all relevant physical quantities in a dimensionless form. To this 
end we introduce the dimensionless photon energy and magnetic field strength 

e = TiLo/nieC? , (1) 
B = B/Bcr = es (2) 
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Figure 1. Geometry for electron-photon scattering. The left (right) panel refers to the situation before (after) scattering; all relevant 
angles are indicated. The direction of the spawned photon with energy e" is also shown (note that this is only allowed for intermediate 
states with quantum number n 2). The figure has only illustrative purposes and here all quantities are referred to an observer at rest 
with the star. 



where ui is the photon angular frequency, B the magnetic field and Bcr = rn^ c^/eh = 4.414 X 10^^ G the critical QED 
field, at which the fundamental cyclotron energy huiB equals the electron rest mass energy. Similarly, we denote with E the 
particle energy expressed in units of meC^ , and with p the component of the particle momentum parallel to the magnetic field, 
expressed in units of meC. 

The process of electron-photon scattering can be schematized as follows: 

(i) an electron in the lowest Landau level at a certain position P in the magnetosphere; the initial spin orientation is 
necessarily down; 

(ii) an incoming photon with energy e which propagates in a direction forming angles 9 and if with respect to the magnetic 
field direction in P (see Figure [T}. The incoming photon may have either linea]|f| polarization s — 1, also indicated by ||, i.e. 
parallel to the magnetic field {ordinary photon) or s — 2, also denoted by _L, i.e. orthogonal to B {extraordinary photon). 

(iii) an excited electron which occupies an intermediate Landau level with quantum number n > and energy E„. 

(iv) a scattered photon with energy e', direction 9', Lp' and polarization s' = 1,2. 

(v) a recoiled electron in the Landau level £ < n, with energy Ee — 1 + e — e' and parallel momentum p'. If ^ > the spin 
of the recoiled electron can be either up {s pin-flip transition) or down {n o spin-flip transition). In the following sections these 
two possibilities will be denoted, following iHarding fc Daughertvl l|l991l ). with / = 1 and f — respectively. 

(vi) If ^ > electron de-excitation to the ground level is accompanied by the emission of one or more photons (Landau- 
Raman scattering). 



In order to derive the releva nt expressions at reson a nce, w e st art from the fully relat i vistic magnetic Compton scat 



Daughertv fc Hardind l|l986l ) and iHarding fc Daughertvl l|l991f ). We do not refer to the 
Gonthier et al.l ( 2000l ). because our final goal is to derive expressions that can be used in 



tering cross section as derived by 
ultra-relativistic limit discussed in 
simulations where different electron population s are present, from non-rela t ivisti c to mildly and ultra-relativistic. The cross 
section derived bv lDaughertv fc Harding l|l98d) and IHarding fc Daughertvl (|l99ll ) includes excitation of the electron to arbi- 
trary Landau levels. However, since we are interested in computing spectra up to a few hundred keVs, scattering occurs (for 
all resonances) where B < 1. For such fields the probability of exciting the second Landau levels is already modest, although 
non- negligible (see ^and Figure [2|, and becomes even smaller (typically by a factor > 10) for n > 2. For this reason we will 
assume that the scattering process occurs only via intermediate states n 2 and, consequently, the final level of the electron 
(i.e. after the scattering) can be either ^ = or ^ = 1. In the latter case the electron remains in an excited state and the 
rapid de-excitation to the ground level is accompanied by the emission of a new photon. The relevance of this process (photon 
spawning) is discussed in Section [5| 

In the following sections we will make use of four different reference frames: a fixed frame (LAB), centered on the neutron 
star; a frame in which the electron has vanishing parallel momentum before the scattering (the electron rest frame, ERF); a 
frame in which the recoiled charge has vanishing parallel momentum (ERF*); a frame in which the electron is at rest in its 
virtual excited state n (ERF**). Unprimed (primed) variables refer to quantities before (after) the scattering. 



^ In general photons propagating in the "vacuum plus cold plasma" are elliptically polarized. However, in our applications, particle 
density is low enough ( < 10^^ cm~'^) to mak e the plasma contribution to the dielectric tensor negligible. The two normal modes are then 
linearly polarized (e.g. (Harding fc Laill2006h . 
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3 RESONANT CROSS SECTIONS FOR PHOTON-ELECTRON SCATTERING 



In strong magnetic fields the transverse momentum of a charge is not conserved upon scattering since the field can absorb 
or add momentum to the particle. The charge momentum parallel to the magnetic field, p, instead, is related to the parallel 
components of the photon momentum before and after the collision by the usual conservation law, which, in the ERF reads 



e cos 6* — e' cos 0' . 



(3) 



After the collision the particle is left in a state with principal quantum number £ ^ and energy E'l — -y/l + p'^ + 2£B. 
Energy and parallel momentum as measured in the ERF and in the frame where the recoiled electron is at rest, ERF*, are 
then related by the Lorentz transformations 

Vl + 2W = Ee = -i'{E'e~0p') and = p* = -f' {p - (5' E'e), (4) 

from which the velocity and energy of the excited particle in the ERF follows 



7 = 



y/l+p'^ + 2eB 



Finally, parallel momentum and energy conservation yields the final photon energy 

sin^ e + 2e- 2£B 



1 



^(1 



i9 cosf) + ^1 + 2e cos6l'(cos6" - cos 6^) + ( cos 6 - cos6l')^ + 2£B sin^ 9' 



(5) 



(6) 



A strong magnetic field greatly afi'ects the scattering cross section, but its effect is substantially different in the case 
the photon electric vector is parallel or perpendicular to B. For frequencies b elow ljb the cross section of the extraor dinary 
mode (_L-mode) is greatly suppressed with respect to the Thomson value, ctt (|Canuto. Lodenguai fc Rudermanlll971f) . Both 
polarization modes manifest, however, a singular behavior at nearly regular frequencies ujn, which correspond to excitations 
of the quantum Landau levels. In the frame where the electron is initially at rest and using the dimensionless variables defined 
in !j2]the relativistic second order cross section for the transition from the ground to an arbitrary state £ is 

2 



do-.. 



Sar e (2 



e - e') exp [-(e^ sin^ 6 + e' sin^ e')/2B] 



n = 



) 



r(2) 



exp (2i<l?) 



(7) 



(see equation [11] of iHarding fc Daughertvlll99ll where some misprints have been corrected). Here di}' = sin 0' d9' dip' , the 
phase <I> depends on the difference ip — (p' , and the complete expressions for the functions F^'']_ can be found in the Appendix 
of Harding & Daughcrty (199ljl. The major complication in equation ([7]l derives from the presence of an infinite sum over 
all intermediate (virtual) states with principal quantum number n. In addition, the F^'']_ are complicated complex functions 
of B, £, e' , 9,8', ip , ipj . They depend also on the initial and final photon polarization mode, s, s' , and on the spin orientation 
of the electron in the intermediate state, i.e. spin-up (labelled by the index +) or spin-down (index — ). Finally, if ^ > 0, the 
form of these functions is different if the orientation of the electron spin in the final state is up (/ — 1) or down (/ = 0). If 
the final state is the ground state only spin down is allowed. 

All the F^^l_ exhibit a divergent behaviour at the resonant energies 

2nB 



1 + \/l + 2nBsm'^e 



(8) 



while the F^l_ , which are obtained from by plane crossing symmetry replacement of variables (see again lHarding fc Daughertv 
remain instead finit^. Since, as mentioned earlier on, we are interested only in resonant scattering, this leads to a 



19911 ) 



(2) 



major simplification. Neglecting all non-resonant contributions in ((T)) when e ~ e„, which amounts to disregard both the F^ 
and the non-resonant terms in the F^^l_, makes it possible to express the cross section as a infinite sum of separated terms 



da 



(n) 



dil' ^ ^ dn' 

n = l 

where each resonant contribution is given by 

^"^ _ 3<7T e' {2 + e- e') exp [-(e^ sin^ 9 + e'' sin^ ( 



da 



n/2B] 



dn' 



IGvr e 



(1) 



(9) 



(10) 



^ Note that, due to a typo, in equations A1-A5 of lHarding & Daughertv! l ll99l[) . all terms Ne,t' A,,A,, must be replaced by N^t,t' A,,A,, 
* Actually also the F^}^ diverge, but th is occur only at photon energies > 2/ sin 6 which is the threshold for pair production and where 
expression jTjl is not valid any more fe.e. lHeroldlll979l) 
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The quantities F^^l_ are independent on ip and have the general form 
V 2Er, J 



^ n,± 

where 



'l + E„\ T^''{e,e',e,e',B; s,s',f) 



En = Vl + cos2 e + 2nB 



(11) 



(12) 



is the energy of the electron in the intermediate state. As discussed in ^ we restrict our treatment to situations in which 
scattering occur only via excitation of intermediate states n ^ 2, i.e. we only consider the first two resonant terms in 
equation (|10|) : the corresponding expression for T"' are given explicitly in the Appendix. 

Actually, the divergence of the f!^X^ which occurs when the denominator 1 + e — En vanishes, reflects an unphysical 
behaviour of the cross section and is merely a consequence of the fact that the expressions presented so far have been computed 
without accounting for the finite lifetime of the electron in the excited Landau levels. In a realistic situation, according to 
the Heisenberg principle, each excited level has an energy indeterminacy which is proportional to the inverse of the lifetime 
of the electron in that state, i.e. to the decay rate. This is standardly accounted for by introducing the natural line widths, 
which amounts to perform the substitution 



En 



En-irl/2, 



(13) 



in the denominator of equation (|lip : here the are related to the relativistic decay rate out of the n-th intermediate state 
and their explicit expressions will be presented in ^ With this substitution, the functions (|ll|l become 



F, 



(1) 



n + En \ _ 

\ 2En ) 1 



— ■ (14) 
+ En+iri/2 ^ ' 

A more useful representation is obtained by replacing the Lorentzian that naturally arises from (I14|) . with a 5-function. By 
taking the limit 

ri/27v 



lim 



5{e -En + 1) 



r'lZo (e-£n + l)2 + (r5/2)2 
we obtain, after a lengthy but straightforward calculation, 



(1) 



-En 



En 



n.l\2 



+ 



'T 



n,l\2 



+ 



y/l + 2nB sin^e 



(15) 



(16) 



We note that the factor E„ (en) / \/ 1 + 2nB sm^ 6 in eq. (|16p arises because of the change in the argument of the (5- function, 
from e — En (e) + ltoe — e„. Ultimately it refiects the fact that the "effective" width of the scattering line profile near resonance 
is TnEn{en)/{2^/l + 2nB sin^ 6) (see lHarding fc Daugher"t^ll99ll ). 

Substituting back equation (|16[l into equation l|10[l and upon a trivial integration on ip' , brings the n-th resonant term 
lfTO)l in the form 



da 



d cos 9 ' 
where 



= v'^i^,{B,e,e')&{e-en). 



(17) 



V 



(n) 



(1 + Enf e'{2 + e- e') exp [-(e^ sin^ 6 + e'^ sin^ ^' 



)/2B] 



16 



+ 



EnV^ + ^nBsin^e e[l + e - e' - (e cos - e' cos 6") cosO'] 

and e is given by equation ((6]|. In the previous expression one can safely put e = en since 2?^^^/ in multipled by S{e — e 
The n-th resonant term in the angle-integrated cross section is 

= S[%(B,e-J,me-en), 
where the integral 

1 



(18) 



(19) 



si%{B,e-,i,f) = 



v[%d{co^e') 



(20) 



can be evaluated numerically. The n-th resonant term of the total cross section is then obtained by summing equation (|19p 
over the polarization states of the scattered photon, s', and over the energy levels I, and the spin states, /, of the recoiled 
electron. 

When n = 1, the terms -S*^!^^, and Ti'f}^^, are not associated with a choice of different quantum states t and /, because 
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-2-10 1 2-2-10 1 2 

log B/B^^ log B/B^^ 

Figure 2. The resonance factors for the first (upper curves, sohd red fines) and second (lower curves, dashed blue lines) intermediate 

(n) 

Landau levels as a function of log-B; the S), are in unit of the Thomson cross section. Different curves are labelled by the value of 
the angle between the incident photon direction and the magnetic field. Left: ordinary photons. Right: extraordinary photons (note the 
weaker angular dependence in the latter case). Clearly, for a photon propagating parallel to the magnetic field only the first resonance 
(2) 

exists and Sg vanishes. 



in this case it necessarily follows that £ — and / — 0. Therefore, for a photon with given initial polarization s, the first 
resonant term of the total cross section is simply given by the sum of two contributions, i.e. 



.(1) 



(21) 



Instead, the second resonant term of the total cross section, corresponding to the intermediate state n — 2, must account for 
both final states £ = 0, 1 and, if £ — 1, for two diff'erent orientations of the electron spin. Its expression is then 



Si%{B,e;£ = OJ = 0) + S^Z2iB,9-£^0,f ^0) 



(2) 



+ 



Si%{B,e;£=l,f ^0) + S^Z,{B,6-£=1J = 0) 



(2) 



(22) 



e;£=l,f ^1) + SlZ2iB, £ = 1, / = 1) 



(2) 



}5(6-e2) = sf\B,9)5{e-e2). 



Finally, the total cross section is obtained by adding the different resonant terms 



as{e,B,e) = ai^\e,B,e) + a^'>{e,B,e) = y ' S';'>{B,e)5{e-e. 



(2)/ 



2 



n = l 



(23) 



The eight functions •S'^'^^z which appear in equations (I21|l - p3| ) are the fundamental ingredients for any computational 
investigation of photon diffusion in strong magnetized plasmas. In particular, when performing Monte Carlo simulations, 
proper combinations of these expressions can be used to evaluate the probability of transition between different states. 

Figure [2] illustrates the dependence of the functions Si^^ and on the magnetic field strength for both polarization 
states and different incident photon direction. The non-relativistic limit tangent to the low _B-field limit of each curve and 
takes the form 

3. 
8q 
37r 



(1) 



^(2 + B cos^ ( 
8a 



-(2 cos^6' -I- B) 



j(2) 



D ■ 2 a 

— B sm , 

2a 



j(2) 



2a 



B sin e{B+ cos^ 



where a = 1/137 is the fine structure constant. In the ultra-relativistic limit no simple analytical expansion can be derived. 
Approximate expressions, which hold for 0/0, are ^i^^' ~ 5''^' oc 1/(1 + 2VB sin^ 6 ) and 5'|'j^' ~ 5''^^' oc 1/(1 + VB sin^ 9). 
As it can be seen, relativistic corrections are already ~ 10% at Z3 ~ 0.1 Bcr for photons with large values of the incident angle 
9 (see Figure [2]). Furthermore, we may note that the number of excitations of the second Landau level, totally negligible in 
the weak B limit, becomes sizable in moderately relativistic regimes. 

This result is more evident in Figure [3] (left panel) , where we show the ratio between the second resonant term and the 
total cross section in the case of unpolarized incident photons, obtained averaging over ordinary and extraordinary modes 
(s = 1, 2). As it can be seen, in strong magnetic fields a non-negligible fraction of collisions can excite electrons to the higher 
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Figure 3. Left: The ratio 5(2)/(5(l) + 5(2)) versus the magnetic field strength for three different scattering angles. The functions 5'"^ 
are obtained by summing over the two initial photon polarization states, i.e. S^"' = 5^"' +52™'. Right: The 2^1 transition probability 
versus magnetic field strength. 



Landau level (n = 2). It is also worth noticing that, when n = 2, de-excitation of the electron to £ = 1 state is generally more 
likely than the direct transition to the ground level, as it can be seen from the right panel in Figure |3l Only for magnetic 
field strengths largely exceeding the critical quantum limit, the two transitions are have comparable probabilities. This has 
important consequences on the transfer problem, since it implies that collisions involving the intermediate state n = 2 lead 
more frequently to the creation of extra photons as the electron returns to the ground level. In other words, the resonant 
magnetic scattering may not conserve the total photon number. 



4 TRANSITION RATES 



For the sake of completeness, in this section we present the explicit expressions for the transition rates and for the natural 
line widths FJ. The differential transition rates (nu mber of decays per second an d per sterad) for electrons with initial spin 
up ( +) or down (—) were derived in a general form by Daughertv &: Ventura ( 19781 ) and Herold. Ruder fc Wunneil ( 1982 ) (see 
also iLatallligiil : iBaring. Gonthier fc Harding |2005| ). Here we will use the derivation presented bv lHerold. Ruder fc Wunneil 
by focussing only on those cases which are of immediate relevance to the present investigation. In this respect, we note 
that, while performing a Monte Carlo simulation, the transitions rates have a twofold role. First, as discussed in §[3l they fix 
the natural width of the Landau levels, and therefore enter directly into the numerical computation of the resonant scattering 
cross section. Second, as we will see in §[5l they determine the properties of the photons emitted by successive de-excitations 
when, in Landau-Raman scattering the electron is left in an excited state £ > 0. In practice, in the particular case we are 
discussing (i.e. only levels up to n = 2 can be excited), it is the differential transition rate dR^'^" /dQ (see equation 1271 below) 
that fixes the properties of the photon which is emitted through the decay of an electron from £ = 1 to the ground level. 

The cross sections discussed in the previous sections have been computed in the frame in which the electron is at rest 
before scattering (ERF, see and the parallel component of the momentum in the intermediate state n is p = e„ cos 9. 
However, in this section we find more convenient to proceed by first introducing the transition rates and line widths FJ 
as computed in the reference frame ERF** in which the electron is momentarily at rest in the virtual state n, and then to 
transform them back in the ERF. In fact, in the ERF** the equations governing the de-excitation and the consequent photon 
emission take their simpler form. Furthermore, in this section we will indicate with E„ and Ei the energy of the electron in 
the state n, i, and with te and 9e the energy and direction (wrt the magnetic field) of the emitted photon, respectively. All 
these quantities are computed in the ERF**. 

In the ERF**, the energy of an electron which is excited to the level n is therefore 



E„ 



Vl + 2nB . 



(24) 



By applying again energy and parallel momentum conservation we obtain the energy of the electron de-excited to the level £ 
and that of the emitted photon 



Ee = En - €e ^ \/l + 2£B + ei cos2 
2(n-£)B 



En + yT+2B (n cos^ 6e + £ sin^ 9e ) 
respectively. 



(25) 
(26) 
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The differential rate for transitions between a generic level n > and the ground state £ ■- 

I 1982ll . takes the form 



amc 

AttU {n - 
2nB 



€e exp {—Z) Z" 



ly.En {En - e.sin^l 

- ee (£n ± 1) 



2nB - (£n =F 1) 



{sin^6»e,0) + 



(COS 



,1) + 2e^ 



COS 



(sm 



,0)} 



: 0, as given bv lHerold. Ruder &: Wunnei 



(27) 



where Z = e\s\v? 9e/2B. The angular factors in "{ )" refer to two possible linear polarization states s = 1, 2. The rate for 
transition from n = 2 to ^ = 1 is 



dRi 



€e exp {—Z) 



dn 



[4B-e, {E„ =F l)](^-2)72 + 



'iTvh En {En - Ee slu^ Oe ) I 

+ 4B - {En ± 1) + 45 {2Z - 3) (sin' e,,0) + 
+ [4B-ee{En T l)]^V2j (cos' 61^,1) + 2el cos 



\4B 



ee {En ± 1)] {Z - ly - 

6»e(^^-2Z + 2) {sin'6»e, 0) | 



(28) 



Since we are considering only the first two resonances, the two previous expressions are all we need in order to evaluate the 
natural widths and photon emission following Landau- Raman scattering (see 

Total rates can be obtained by integrating over the solid angle and summing over the two polarization states s of the 
emitted photon, and in general take the form 



= 27r 



2 



d(cos 6^ 



■ d(cos 6s 



(29) 



The reciprocal of the transition rate gives the mean lifetime of the electron in the excited Landau level. Then, due the 
Heisemberg principle, each term i?"^* is proportional to the energy width for the particular transition n ^ £. The total 
energy width of the n-th level is therefore obtained by summing over all contributions associated with level n 



h 



(30) 



As previously stated, these expressions are valid in the reference frame in which the excited electron is at rest, ERF**. 
In order to obtain the energy widths that must be used in the expressions given in § [S] we need to perform a Lorentz 
transformation back to the ERF. This gives 



-I + 



VI + 2nB 



Vl + En COs2 



-pn 
J- ± ) 



(31) 



* + 2nB 

where and 6 are the incident photon energy (equation [8]) and direction. Figure |4] shows the dependence on the magnetic 
field strength of the total energy widths for the first two Landau levels, F^^ , in the range 10"=* B/Bcr ^ 10^ For weak 
field s they are proportional to -B^ (in particular Fj_ ~ 4ai3^/3), except for the spin fiip transition for which F+ ~ 2Qf_B"^/3 



(e.g. iHerold. Ruder fc Wunnei 1 19821 : iLatal 1 19861 ). Similarly to what happens for the cross sections, remarkable deviations 



from the simple, non-relativistic power-law dependence appears at B > 0.1. We note that in the st rong field lim it all energy 
widths have a similar dependence on the field strength, ~ 0.4q:V'B, a result similar to that found by Latall | l986l) for the total 
transition rate. The limiting behaviours of F± can be usefully checked against those derived by Baring. Gonthier fc Hardin j 
( 2005h who obtained an analytical form (in terms of a series of elementary functions) fo r the spin-dependent transition rate s 
R"-^°. In fact, for n = 1, it follows from eq. (|30} that r± = (ft/meC^ )-R+"° . As ri oticed bv lBaring. Gonthier fc Hardinj (|2005l '), 
their weak field limit coincides with those derived (analitycally) by Latall ( 19861 ) and discussed above. In the ultra-relativistic 
limit, on the other hand, one gets r_ ~ r_|_ ~ (1 - l/e)(B /2)^/^ ^ 0.45a/B, close to t he B > 1 limit given above. A 
direct numerical comparison of our result for Fj- with that of Baring. Gonthier fc Hardin j (2005) shows that the fractional 
difference is below 10% for B < 100, with the largest deviations appearing at large B. 



5 PHOTON SPAWNING 

As discussed in the previous sections, some scatterings occur at the second resonance and excite the electron in the intermediate 
state n = 2. In order to have a significant number of second harmonic excitations the magnetic field must be sufficiently intense, 
say B > 0.05 {B > 2.2 x 10^^ G; see Figure |3]). Moreover, the incoming photon must have an energy e = e2 « 2B. The two 
conditions together require that a significant number of photons with energy hu > 50 keV are present in the magnetosphere. 
We do not expect such high-energy photons to be emitted directly by the stellar surface, which temperature, as inferred by 
observed X-ray spectra, is ~ 0.5-1 keV, but it must be noted that, if electrons are relativistic, the photon energy as seen by 
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-1 h H -1 




-5-2-10 1 2-3-2-10 1 2 

log B/B^^ log B/B^^ 



Figure 4. Energy widths of the Landau levels (in units of m^c^) for the two spin orientations as a function of the magnetic field strength. 
Left: n = I. Right: ra = 2. In each panel the full lines refer to the widths as computed in the excited electron frame [oq. Il30ll ] and the 
dashed lines to the same quantities but referred to the ERF [eq. | |31| |] for cos 6 = ±1, i.e. in the case in which the incoming photon is 
parallel to B. 



the moving particle will be amplified by a factor ~ 7. Besides, repeated collisions populate the hard tail of the spectrum and 
may provide photons of the required energy. 

Once an electron is excited to the second Landau level, it is more likely left after scattering in the £ — 1 rather than 
in the £ — state (see § |3]and Figure [3]). This implies that an additional photon will be produced by the prompt radiative 
de-excitation 1 ^ 0. In general, the transitio n rate depends also on the rate at w hich collisions populate the £ — 1 level. The 
excitation rate due to Coulomb collisions is ( Bonazzola. Hevvaerts fc Pueet 1979[ ) 



Rc 



3 X 10' 



Vl02icm-3y 



T 



-3/2 



1.6 X 10' 



(lO^icm-s) ^ 



3/2 



where the second (approximate) eq uality follows from the re quirement that it has to be kT ~ huB for the collision to excite 
the first Landau level. As noticed by Arava Hardin3 1 1999I ). in a strongly magnetized, low-density plasma Rcoii is negligible 
when compared to the radiative cyclotron rate 



Brad ~ 10'* B" s-i (32) 

where the exponent q varies in the range 1/2-2 (see |4|. The population of the excited Landau levels is, therefore, solely 
regulated in the present case by magnetic Compton scattering. The large value of Rrad justifies the assumption that electrons 
remain in the ground level. In fact, the typical time between two scatterings is « (riecrc)"^ « L/{tc) ~ 10~^ s ^ 1/Rrad, 
where the lengthscale L is a few star radii and the scattering depth r ~ 1 (see paper I). In a sense, one may view Raman 
scattering at the second resonance, and the ensuing, quick radiative de-excitation to the £ — level (photon spawning) as 
more akin to double Compton scattering (e.g. Gould 1984) rather than to absorption followed by the emission of two photons. 

As discussed in § 3] the de-excitation of an electron which is initially at rest in the £ — 1 level with energy E = a/I + 2-B 
produces a photon with energy 



E - ^/E^ ~ 2B sin" ^ 2B 

sin^^e VI + 2B -h VI + 2B cos" ' 

as computed in its rest frame. The rate at which photons are emitted is given by equation (|29|) with n = 1, £ = 0, and, when 
performing a Monte Carlo simulation, the differential form (|27|) can be used to derive, on a probabilistic ground, the direction 
of the emitted photon (i.e. the angle Q^, emission is isotropic in 0^). Once again, we stress that care must be payed to the 
frame in which these quantities are evaluated. In fact, the emission angle and the photon energy introduced above, are 
both referred to the rest frame of the recoiled electron, while in applications they need to be evaluated in the stellar frame 
(LAB). This is readily done by means of a Lorentz transformation. In this respect, we note that the parallel momentum of 
the recoiled electron is in the LAB 

V ~ 7/3 e /i — e' (34) 

where e and /i are the photon energy and the cosine of the propagation angle with respect to the magnetic field, both 
measured in the stellar frame, and, again, unprimed (primed) variables refer to quantities before (after) the scattering. The 
corresponding velocity and the Lorentz factor are /?' = p7^1 +p'" -I- 2B and 7' = ^/TTp^^~^2B / VTT2B respectively 
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log B/B^^ 
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0, 

o 

0, 
0, 
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o -0. 





0.0 
cos 9' 



Figure 5. Le/i; the ratio (Taha/'^s versus tlie magnetic field strength for ordinary (upper panel) and extraordinary (lower panel) photons, 
and three different directions of the incident photon, cos9 = 0,0.5,1. Right: the angular distribution of scattered photons (upper 
panel) for a given incident photon direction 6 = 30°. The lower panel shows the fractional difference between emission and scattering, 
{Qe — Qs)/Qs - Each curve is labelled by the magnetic field strength. 



(see from which we obtain 

ML) _ cos Qe + P' 



(35) 



1+13' cos Oe 

and 

ei^' = 7'ee (l + /3'cosee) ; 

here the index (L) denotes the photon energy and propagation angle as measured in the stellar frame to distinguish them 
from the same quantities but referred to the recoiled electron frame. 



(36) 



6 ABSORPTION-EMISSION VS. SCATTERING 

As noted by Harding fc Daughertv 1 1991 ) and Arava fc Harding ( 19991 ). in the non-relativistic regime the on- resonance second 



order scattering process y + e —> ^' + e' is well approximated by the sequence of two separated first order processes: photon 
absorption, 7 + e — > e*, immediately followed by photon emission e* —> e' + 7'. Because of its much simpler form, one would 
like to use the cyclotron absorption cross section instead of the cumbersome QED scattering cross section in numerical work. 
It is therefore of interest to explore the conditions under which absorption plus emission reproduces relativistic resonant 
scattering with reasonable accuracy. 



19781 ) 



In the ERF, summed over the final spin states, the angle-averaged absorption cross section can be written as IjDaughertv fc Ventura 



(II, 



Snar exp(-Z)Z"" 



4" {n-iy.y/l + 2nB sin^ 



nB 



,0) 



- En), 



(37) 

el sm^9/2B, and 



where the angular factor { sin^ 6, 0) refers to the two linear polarizations (||, _L) of the incident photon, Z 
the cyclotron harmonic energy e„ is equal to the scattering resonant energy (see equation [8]). 

The reliability and the limits of validity of this approach can be assessed by directly comparing equation (|37|l with the 
scattering cross sections (equations [21] and [23]), since these expressions are all explicitly written in terms of a 5- function 
of the same argument. Figure [5] (left panel) shows the ratio of the absorption to the scattering cross sections for n = 1 and 
both ordinary and extraordinary photons in the range —2 ^ log B ^ 2 and three different directions of the incident photon 
(corresponding to cos 9 = 0, 0.5, 1). The agreement is very good up to B ~ 0.1 for both polarization modes and deviations are 
within ~ 20% up to B ~ 1. However, this does not imply that the two descriptions of the photon-electron interaction are to 
be regarded as equivalent for B < 1. To claim this, one should prove that they produce also the same angular redistribution 
of radiation in the final state. The angular distribution of the scattered photons is given (in the ERF) by equation (|18p once 
B and 9 axe fixed, while that of emitted photons by equation (|27|l . The latter is evaluated in the zero-momentum frame of the 
excited electron, exhibits forward-backward symmetry, and lacks any information on the incoming photon. However, in the 
case under examination, the electron has been excited by the absorption of a photon with parallel momentum e„ cos 6. In order 
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to co mpare the two angular distributions the transition rate must be properly transformed to the ERF IjArava fc Hardine 



1 19991) . or, alternatively, it should be computed ab initio in the ERF, by exploiting the general expressions derived in iLatal 
l|l986l ) and valid for an arbitrary value of the electron parallel momentum. Since p — €„ cos 6, this clearly introduces a 
dependence on the incoming photon kinematical quantities in the transition rate. 

Figure [5] (right) shows the angular distribution Qa{9') of the scattered photons (upper panel) for n = 1, a representative 
value of the incident photon angle, 9 = 30°, and different values of the magnetic field, 0.01 ^ _B ^ 10. Qs has been computed 
by normalizing equation (|18p to its maximum value, which occurs at 0' = and is the same for any B. Photons are more and 
more forward-scattered as the field increases. We verified that, as expected, this asymmetry strongly increases for 6 ~ and 
6 ~ 180°, while the curves become symmetric for 6 = 90°. The lower panel shows instead (again, for n = 1) the fractional 
difference between Qs{6') and the analogous quantity for emission Qe{0') (here 9' = 9^), obtained normalizing the Lorentz- 
transformed transition rate (|27p . It is interesting to note that in this case the forward beaming is entirely due the relativistic 
Lorentz boosting in going from the ERF** to the ERF. As the plot shows, there are significant differences in the angular 
distributions even at magnetic field below O.lScr. This means that photon-electron scattering can be safely treated as the 
superposition of absorption and emission only for resonant photon energies < 10 keV. 

Similar results are obtained for transitions involving higher Landau levels, although in this case sensible differences 
arise even for lower magnetic field strengths. However, when n > 1 major complications arise because of the necessity to 
discriminate transitions toward final levels different from the ground state, and to account for the possible changes of the 
photon polarization and electron spin orientation. 



7 SCATTERING PROBABILITY 

Let us assume that a photon propagates in a strongly magnetized medium populated by electrons with a one-dimensional 
relativistic velocity distribution along the magnetic field direction 



(38) 



where fe ~ n~^dne/d{'yf3) is the momentum distribution function and ne(r) is the charge number density. As discussed in 
^ all electrons can be taken to be initially in the ground Landau level. Having in mind the results obtained in Section [3l 
equation (|23|) . the optical depth after a photon has travelled an infinitesimal distance di is 



dTs 



d£ 



2 

n = l 



df3n,ir) 7^ /e(f, 7/3) (1 ~ /J^) S: S [7(1 - /3/i)e(^> 



(39) 



/9„ 



in which [l3min,Pmax] is the charge velocity spread and jj. is the cosine of the propagation angle with respect to the magnetic 
field in the stellar frame (LAB). The latter quantity is related to the same angle as measured in the ERF by the usual 
transformations 



cos U = 



l-/3/i 



7(1 - /3/i) 



(40) 



The factor (1 — f3fj,) in the integral (|39|l appears because of the change of reference between ERF and LAB, and, for the sake 
of clarity, the argument of the 5-function has been explicitly written in terms of the energy e'^' measured by an observer at 
rest in the LAB. 

The integral (|39|l can be readily performed exploiting the well-known properties of the (5-function. Denoting by /3^"^ {k — 

' e„ = 0, the (5-function in energy can be transformed 

(41) 



1, 2) the two roots (for each n) of the quadratic equation 7(1 — /3/x)e'^' 
into a (5-function in velocity by 

' <5(/3-/3i")) 



5[7(1 - /3M)e<^' 

where 
9(0) = 



E 



v/7^(l-/3M)' + 2nB(l-^2) 



(n) 



1 



(42) 
(43) 
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nB 

En 



el (1 - M^) 
2nB 



(44) 



It is worth noticing that with this transformation the role played by the resonant photon is replaced by that of a pair of 
resonant electrons selected among all charges of the distribution. Accordingly, the elementary optical depth (|39p for a photon 
with initial polarization state s = (1,2), becomes 



dTs 



Uedl 



(")\ q(") 



sr(B,e(M,/3<"')) 



n=l = l 



Iff I 



(45) 



the quantities As"' are implicitly defined by equation (|45p , and we made evident the dependence of the resonant factors Ss"' , 
which are computed in the ERF, on the LAB variables. 

Once energy, polarization and initial photon direction are fixed, equation (|45|l can be integrated numerically along the 
optical path until a scattering, if any, occurs. As it is apparent from equation (|42|) . scattering may occur only when the roots 



/3^"' are real, i.e. only where h(r,fi) 



2 I 2 

; m + fJ. 



1^0. Since the function h depends only on position and photon direction. 



at every point in the magnetosphere the previous condition discriminates those pairs of energy and angle for which scattering 
is possible. 



8 CONCLUSIONS 

Recent models of spectral formation in magnetars called renewed attention on electron-photon scatte ring in the presence of 



ultra -strong magnetic fields. The complete expression for the QED cross section is known since long (e.g. lHarding fc Daughertv 



19911 ) but its practical application in the relativistic regime is numerically challenging. In many astrophysical problems, 
including the one which motivated us, it is reasonable to assume that scattering occurs only at resonance, i.e. when the incident 
photon frequency equals the cyclotron frequency (or one of its harmonics). Restricting to resonant scattering introduces a 
major simplification, and here we presented explicit expressions for the magnetic Compton cross section in this particular 
case. Our main goal has been to provide a complete, workable set of formulae which can be used in Monte Carlo simulations 
of photon scattering in strongly magnetized media. 

Our results are fairly general and can be applied to resonant photon scattering under a variety of conditions. In particular, 
no assumption is made on the field strength. Having in mind applications to spectral modelling in the ~ 0.1-200 keV range, 
resonant scattering necessarily occurs where B/Bcr < 1. Under these conditions, and although the expressions we derived are 
still valid for B/Bcr > 1, we restricted to the case in which the electron is excited at most up to the second Landau level. 
We find that deviations from the non-relativistic limit in both the first and second resonant contributions to the cross section 
become significant for B/Bcr ^ 0.1. The probability that scattering occurs at the second resonance, which is negligible below 
B/Bcr ^ 0.01, becomes sizeable at higher B and, depending on the scattering angle, it can be up to ~ 30% for B ~ Bcr- In 
case the second Landau level is excited, its is more likely that the recoiled electron is left in the first than in the ground level 
with the ensuing emission of a new photon (spawning). Using our results for the cross section together with known expressions 
for the transition rates, we checked under which conditions resonant Compton scattering can be treated as the combination 
of two first-order processes, photon absorption followed by emission. While the scattering and absorption cross sections differ 
by at most ~ 20% for B < Bcr, the angular distribution of the scattered/emitted photons shows deviations already at 
B ~ O.lScr. Finally, having in mind the implementation in a Monte Carlo code (Nobili, Turolla & Zane, in preparation), we 
presented an explicit derivation of the scattering optical depth along the photon path. 
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APPENDIX A: 



The expressions below explicitly give the functions T"^^{e, e' ,6,6' , B; s, s' , f) introduced in equation (lllf) for the relevant 
values of n and £. In the following, to simplify notation, only the dependence on / is shown in the argument: / = refers to 
the no-flip and / = 1 to the flip case. 



Ti-*°(0) 



(2 + 6)2 

1 

(2+^ 



G 

— (Xi COS6I-X3 sine*) + sin6'' (tVIo sin6'-7\/(2 cose*) 



He' sin6>' 
G 



2B (2 + e)2 
ee' sinS sin 6'' 



(2B7M2 +Moe^ sine cosd) + {2BMi +Mae^ sin 6 cos f 



(A^i cos6' — A^3 sin6') + sin6'' (yVfo sm9 — M2 cos( 



He' sin 6*' 



2B (2 + e)2 [ 
^e^e'sine (G 



{ABM2 ^ Mae sine cose) + {'iBMi +M:i£ sin6 cost 



453/2(2 + e)2 

"2B(2 + e-2e') 



|^(2B-e'^ sin^e')(Xi cos^-Xg sin6l) 



+ smf 



{Mo sin 6* - M2 cost 



\/2" e sin 6* ( H e' sin 61' (4B - e'^ sin^ e' ) 



8B3/2(2 + e)2 1 2BD. 
+ {4B-2e'^ sin^ 6' ) {4BM 1 + X 3 sin 6* cos 6^) + 



{ABM2 + e'^Mo sin6 cos 6) + 
sin^6'{2 + e) 



n 



(ABMb + Mi e sin 6 cos t 



i^He sine' {M2 cose-TWo sine*) + (2 + e)(2B - e'^ sin^ e')(7Wi cos^-TWg sine*) + 



25^(2 + 6)2 

+ ne'^ sin e' {M5 cos6-Mi sin 6)^ 

e e'^ sin ^ sin 6' r ^ , . „, , , ^ , ^ . , o ^ . 
- gi^g2 (2 + e)2 l*^ ™ ^ (4B AI5 + X4 e cos e sin 6l) + 

+ [AB-e sin^6l'(2 + e)] (4BX2 + Xo sine cos6l) j 



where 

G = (2 + e) cos e' — e cos 
and 

, sin £» sin 6' 

Mo = 
X2 = 
7W4 = 



// =e(4 + 2e-e') cos6l-e'(2 + e) cose*' , Q^2 + e-e 





- cos 6 sin e' 
sin e' 

sin^ cos 9' sinf 




Ml = 
M3 = 
Ms = 



cosy cosy cost 
- COS 61' 1 

— sin 6* COS 61' 




sm f 




COS ( 

-1 



In each expression the functional form for given values of the photon polarization, s, s' , is recovered inserting the matrix 
element Mk {s, s') . 



